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DECUS Program Library Write-up 


DECUS NO. 8-566 


Siy 

SUB 8IC-# Jan. 22, 1972 

(1) SAME ; PARTL 

Purpose : To evaluate the partial fraction 

expansion o. a rational function 
N(s)/t)(s). that has real coeffi¬ 
cients and D(s) are written in 
linear of quadratic factors* 


(2) CALLING SEQUENCE : 

CALL PARTL(N,IPOLE,SP,P,Q,RI 5 AI,BI) 
where: 

Input P&t&t 

* deg D(s) 


N -_„ 

I POLE a 1, means linear factor in D(e) 

a 0, means quadratic factor in D(s) 
A(I), Ia1,2, «.», N coeff* of N(a) 

B(I),I *1,2, .... N+1 coeff* of l)(s) 

SP a (see (4), or (7), or (3) ) 

P 

Q 


( 

( 


If 

ft 


ft 

ft 


f! 

If 


IS 

ft 


M 

It 


It 

!! 


Return Data : 

HI * < 11 

AI « ( ” 

BI * (" 


it 

ti 

ti 


tt tt 

it t» 

it it 


it tt ) 

H ti ) 

II It ) 


(3) ERROR RETURN : 

-Ho build in error check, but if computer over 
flow /underflow the following message will be 
printed on the Teletype printer,, 


"FPNT" ERROR AT LOC 031^2 

(4) SPECIAL CONSIDERATION : 

User should make sure that the following strict 
requirements are followed: 

1. degD(s) a N, at least 2 

2. degN(s), is at most N-1 

3* D(a) must be factored into linear and/ 
quadratic factors of the form: 

p 

a + sp, tr + ps + q 

where sp„ p and q are real nos* 











(5) SUBPROGRAM CALLjED ; 

-none- 

(6) STORAGE REQUIRED : 


pages (octal) 

(7) ALGORITHM & REFERENCES : 
Given ; G(a)» N(s)/D$s) s 


a,**" 1 ♦ 


A. 


n^'J- j-. 1 * 

°1® ^ e 2 s * 


<• © 


a*1 


where D(s) » T|~ (s+1,) Tf (a 2 ♦ Pd a * q .) 

£■1 j= i J 


dor some reals 1^, , q ^ 

and integer (poa) m 


Problem : find 
G(s) - 

Solution : 


a K. 

i rr n 


n-*a 



r? 


A.s ♦ B. 
J d 


>"1* + Pj® ♦ 


For each i, let 1* * sp, and for each 4 
Pj » <3lj be * p and q«. 

Assume the restriction in (4) hold, then 
Algorithms 

( CASE I ): For each i «1,2, ..., m 

DCs) 9 e^s 11 ♦ ♦ e n+ ^ 

* (s ♦ l i )(bjS n ~ 1 ♦ 0 . 06 + b Q ) 


stepl ; Generate sequence of nos» b 1t b 2 
and store away, from 


5 « » o 


> b 


b i * *i 


\ “ ®k " V ^-1 for k »2<3, .... n 
step2 : Construct tw sequences, c^,c 2 , <»•<,, c Q 

and , d 2 , -o,, d n and store away, from 


2 



















Cj a a 1 


1 \ 

o « o o O O W 


c k ° ~ c k-1 1 i +a k 


\ * b k 


for k a2 ? 3» *••, n 
etep3 $ Then calculate and print, 

*i * V d n 

atep4 s If 1 » m. go to Case II, otherwise repeat 
step 1 thru step 4« 


( CASE II ) : For each >1,2, n-m 

D(s) a e^ 14 ♦ *.«* ♦ e n+1 

" <** ♦ S 3 S * < b 1 - ••• + V) > 


stepl : Generate seq« bjjbg, <*•♦» b n _.j and store away 
from 

bj ■ e^ , b 2 a °2 * 

\ * ®k “ p J b k~1 “ q jhk-2 for k*3,4».« ,n~1 

set b n a 0 (Important I!) 

step2 : Construct 2 sequences c^c^, c Q and 

dj ) dg , d Q and store away 5 from 

C 1 * a i d l * b l 

e 2 * * ®2 d 2 ■ - d 1* > J * b 2 


oeoeove** o * 

c k * ■ e k-1 p 3 -c k-2 q J + ®k 


\ " -VlPj * d k-2 q ;| + \ 

for k*3,4,c, n 

atep3 t The evaluate, and print 

J Vn^Ci^ 

rj d n c n-1 “ ‘Sx-I 011 

B. > - - -- 

J Vn-2 “ Vi 
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fltepif : If j a n-m, FINISHED. Otherwise 
repeat step 1 thru step 


References: 

BoO. Watkins, "A Partial Fraction Algorithm," 
IEEE Trans on Automatic Control vol. Oct. 1971, 
pp. **89-491. 


(8) LISTING : 

(See attach) 

(9) SAMPLE : 

(See attach) 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


**** PARTIAL FrACTIO iV EXPANSION 

**** NCS)/DtS> = C At l )*S** tfv-I ) ♦ ... + ACM) )/ 

tEtl)*S**N * ......... + KCN-+1) 

tv IS ATMOST ID 
JAN. 22# 1792 

** REQUIREMENTS t 

1. BEutDCS)) » N* AT LEAST 2 

2. DECattVtS ) ) ® AT MOST tN-l ) 

3* Bt S ) MUST HAVE SIMPLE POLES 
A. FACTORS OF DtS) MUST 3E OF THE FORM: 

S + SP* S**2 + p*S + O 
WHERE! SP* P ft <5 ARE REAL EOS • 

IPOLE * 1# SIMPLE POLES 

= 2# QUADRATIC FACTORS IN DCS) 
SUBROUTINE PARTLtN* lPOLE#SP*P#Q#Hl #AI >31 ) 
COMMON A#E 


DIMENSION A C10)#EC11>*Btl0>*CCll)*Dtl1> 
IFtlPOLE-1) 10*10*20 
10 B(l)eECl) 


DO 30 K«*2#N 

30 3tK)»EtK)-SP*3CK-l) 

CCl)®Atl) 

DtI)®nti) 

DO 40 K®2#N 

D< K)= -DtK-l)*SP ♦ BtK) 

43 CtK>« -C<K»1)*SP ♦ At K) 

Rl=CtK)/DtN> 

RETURN- 

20 3tl)»Efl) 

Bt2)«EC2)-P*BCl) 

31N ) =3 « 

M*N-l 


DO 50 K»3*N1 

53 BfK)=ECK)-P*3CK-l)-Q*3tK-2) 

Cfl)®Atl) 

Ct2)®-CCI)*?+AC2) 

Dfl)®3fl) 

Dt2)s-Dtl)*P+Bt2) 

DO 63 K«3#N 

DtK)e-DCK-T)*P-DCK-2)*Q+B<K) 

63 CCK)»-CCK-l)*P-CtK-2)*3+AtK> 

SS®DtN)*DtN-2)~DtN-l )**P. 
AI®tDtN-2)*CtN)-DtN-l)*CtN-l)>/SS 
BI®tDtNi)*CCN-l )-DCN-l >*CCN))/SS 
RETURN 
END 


) 
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c *** * 

c 

c 

c 

c 

c 

c 


iff 

5 


6 


15 

3ff 

25 

C 


43 

45 

53 

55 

C 


63 


< 5 &MFl£- 


SAXHLK <' 1 
CALL PARTL 

N = DEGREE OK DCS) AfXCST 13 

NSR a n SIMPLE HULLS 

NOD « * QUADRATIC FACTORS IN DCS) 

COMMON A#3 

DlXKtvSlU.N) AC13)#3C11) 

READC2#10 ) N#NSP#N3D 
FURXAT C 315) 

DO 5 1*1 # N 
HEADC2#25> AC I) 

N’laN + l 
DO 6 I *1 # Nil 
hKADC2»25) HCI) 

1P0LF. *l 
I *1 

IF CI-NSP) 30#30#40 
RKADCP#«S) SP 
FORMATCEl 2 *5) 

CALL HArTLCN* 1P0LE#SP#P#0#hl#AI>31) 

WRITECP#25) HI 

1=1 + 1 

GO TO 15 

IPULF. *2 

1=1 

IF Cl-NOD) 50#50# 63 
READC2#55) Pi f) 

FO! AAi C2K12.5) 


CALL PAR TL < u * IPOLL # SP #P * 0 *: 1» A l » JI) 
will TEC2#55) AI #01 
1 = 1+1 
GO TO 45 
STOP 
END 
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2 


4 

b • 

9. 

3 • 
- 2 . 
1. 
5. 
13. 
13. 
4. 
1. 
2 . 
2 . 



1 


2 . 



gS* 4-?S*> 3 i -2 
S*+5S i j-IOS % i\OS+l 



! Mwj 

S+J SIJS-t-2 


-.lPi3 , iv i !E+ni 

.60333E*fll 

• 3003C1E+CJPI -.53000E+31 


K^-l 

G* 

A -cp B s ~X 


4 0 2 

9. 


32*15 

241*31 

125.5 

1. 


mi 

6.33 

59.34 

151 .5 

631 .85 

3* 

25. 

DCs) 

3.33 

25.25 




' 

.50300E«-O1 

•20000E+01 


•40000E+91 

.30003E+01 



f %*+3**iSS +#b*/S -HXX 

ta.es Xteutf 


fj>± 2 f ^gfi 

s 2 +3 *V &03S4aC*f 











1* 

1. 

«*O330030002 

2 . 

2 • 3000000 ( 10 ! 


"FPNT" EUR Ok AT LOC 03152 
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